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ABSTRACT 


Solar  flux,  and  "total*  particle  and  02  density  distribution,  in 
the  upper  mesosphere  and  lower  thermosphere  (70-  to  150-km  region) 
are  described.  Respective  distributions  were  determined  from  solar 
flux  intensities  measured  by  ion  chambers  in  three  rocket  vehicles 
launched  on  21  February  1966  from  North  America  (two  at  Eglin  Gulf 
Test  Range,  Florida,  and  one  at  Fort  Churchill,  Manitoba).  Method 
of  determining  density  distributions  is  given. 

Solar  flux  outside  and  inside  the  absorbing  atmosphere  was  meas¬ 
ured  for  the  following  wave  bands:  1-20A  (X  ray),  1140-1260^  (mainly 
hydrogen  Lyman-a  line  at  1215. 6^),  and  1325-1475)1  ( Schumann -Runge  con¬ 
tinuum).  Distributions  of  02  and  "total"  particle  density  are  based 
on  the  solar  flux  absorption  in  the  sensible  atmosphere. 

Results  of  solar  flux  measurements  compare,  in  general,  with 
results  obtained  elsewhere.  Estimated  errors  vary  from  10  to  100 
percent.  Density  results  also  agree  generally  with  those  of  other 
experimenters  (using  similar  or  other  techniques),  except  for  one 
striking  difference:  a  step  of  02  density  exists  at  113  i  7  km  alti¬ 
tudes.  This  step,  which  to  our  knowledge  has  not  been  reported  else¬ 
where,  is  described  theoretically.  Its  existence  is  also  confirmed 
by  electron  density  data  derived  by  ionosonde  at  each  launch  site. 
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I.  INTRODUCTION 


Experiments  described  in  this  paper  were  undertaken  as  a  con¬ 
tribution  to  an  on-going  scientific  effort  to  study  the  upper  at¬ 
mosphere  in  the  70-  to  150-km  altitude  region  (upper  mesosphere  and 
lower  thermosphere). 

Numerous  factors  make  this  part  of  the  upper  atmosphere  an  in¬ 
teresting  one  to  study.  For  example,  most  modes  of  radio  communi¬ 
cation  are  sensitive  to  density  distributions,  temperature,  and 
motions  (winds  and  turbulences)  found  here.  Also,  density  distri¬ 
butions  and  temperature  in  the  upper  mesosphere  control  the  ven¬ 
tilation  upward  of  the  stratosphere  and  mesosphere--signif icant 
regions  in  the  light  of  predicted  pollution  by  exhaust  gases  of 
supersonic  transport  vehicles  and  rockets.  In  the  lower  thermos¬ 
phere,  density  distribution  at  the  higher  altitudes  fluctuates 
widely  as  compared  with  that  at  the  lower  altitudes  where  density 
fluctuations  are  smaller.  The  strong  fluctuations  at  these  higher 
altitudes  can  strongly  affect  the  temporal  changes  in  orbital  charac¬ 
teristics  of  earth  satellites  whose  perigees  may  be  low. 

Study  of  this  region,  in  which  geophysical  processes  are  many 
and  varied,  has  been  handicapped  in  the  past  because  various  experi¬ 
mental  techniques  that  have  been  effective  in  studying  atmospheric 
characteristics  above  and  below  the  region  proved  ineffective  with¬ 
in  the  region. 

Experimental  data  were  obtained  by  three  rocket  vehicles  (Gro- 
becker,  1967a),  launched  on  21  February  1966  by  the  Naval  Research 
Laboratory  (NRL).  Ttfo  lifted  off  at  Eglin  Gulf  Test  Range,  Florida, 
and  one  at  Fort  Churchill,  Manitoba.  Launch  locations,  times,  set¬ 
tings,  sun  positions,  and  apogee  information  are  given  in  Table  1. 
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TABLE  1.  GEOMETRY  AND  TIMES  OF  ROCKET  VEHICLE  EXPERIMENTS 
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Each  rocket  vehicle  contained  four  pairs  of  ion  chambers:  two 
pairs  sensitive  to  vacuum  ultraviolet  rays  (Grobecker,  1967b)  and 
two  sensitive  to  X-radiation (Grobecker,  1967c).  Molecular  oxygen 
density  and  "total"  particle  density  in  the  upper  mesosphere  and  lower 
thermosphere  were  determined  from  the  chambers’  measurements  of  the 
sensible  atmosphere’s  absorption  of  solar  radiation  along  a  near- 
vertical  trajectory  through  the  region. 

Solar  flux  absorption  inside  the  sensible  atmosphere  and  solar 
flux  intensities  outside  the  sensible  atmosphere  were  determined  for 
three  wave  bands  yielding  the  most  useful  information:  1-20A  (X  ray), 
1140-1260A  [mainly  hydrogen  Lyman-er  (H  Ly-a)]  line  at  1215. 6A,  and 
1325-1475A  [Schumann -Runge  (SR)]  continuum. 

Absorption  of  solar  flux  by  the  0^  molecule,  at  the  measurement 
altitudes  of  the  experiments,  was  determined  at  the  H  Ly-o  line  and 
in  the  SR  continuum.  Absorption  by  oxygen  and  nitrogen  atoms,  either 
free  or  bound  in  molecular  form,  was  determined  in  the  1-20A  wave  band. 

Solar  flux  in  each  wave  band  was  measured  by  two  ion  chambers 
(detectors  A  &  B)  in  each  rocket  vehicle.  Originally,  it  was  planned 
to  determine  solar  flux  in  the  44-60A  wave  band,  and  two  detectors  in 
each  rocket  vehicle  measured  intensities  in  this  band  during  the  ex¬ 
periments.  However,  results  from  these  detectors  were  considered  not 
useful  since  the  44-60A  determination  proved  sensitive  to  errors  in 
1-20A  flux  determinations.  Consequently,  these  particular  results 
are  not  discussed  in  the  paper. 

Throughout  this  paper,  detectors  for  each  of  the  three  wave 
bands  providing  usable  d^ita  are  referred  to  either  by  wave  band  (i.e. , 
1-20A,  H  Ly-a,  and  SR  detectors  corresponding  to  the  respective  wave 
band  ranges)  or  by  the  channel  through  which  measurement  data  were 
telemetered  (i.e.,  Channel  14  detector  for  1-20A,  Channel  15  for 
H  Ly-a,  and  Channel  18  for  SR). 

Vertical  profiles  of  electron  density  were  determined  continu¬ 
ously  by  ionosonde  at  each  launch  site  during  launch  periods. 
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Solar  flux  determinations  are  discussed  in  Section  II  of  this 
paper.  They  are  described  in  terms  of  flux  absorption  inside  the 
sensible  atmosphere  and  flux  intensities  outside  it,  flux  computa¬ 
tions,  flux  as  a  function  of  altitude,  and  error  contributions  in 
the  flux  measurements. 

"Total”  particle  and  02  density  distribution,  based  on  the  solar 
flux  determination,  is  discussed  in  Section  III.  This  is  presented 
in  terms  of  analyzing  solar  flux  absorption  data,  computing  density 
values  as  a  function  of  altitude,  and  interpreting  the  experimental 
data  from  the  geophysical  standpoint. 

Information  related  to  observations  of  solar  emission  constancy 
during  the  experimental  period  is  presented  in  Appendix  A.  Material 
amplifying  density  computations  in  Section  III  is  given  in  Appendix  B. 
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II.  SOLAR  FLUX  DETERMINATION 


A.  FLUX  ABSORPTION  WITHIN  THE  ATMOSPHERE 

Individual  results  of  the  three  rocketbome  experiments  to  de¬ 
termine  the  absorption  of  solar  flux  in  each  of  the  three  wave  bands 
yielding  useful  information  are  given  in  Figs.  1,  2,  and  3.  Flux 
values  within  the  sensible  absorbing  atmosphere  are  shown  as  functions 
of  altitude,  determined  from  the  ion  current  of  each  of  two  detectors 
(A  &  B)  turned  toward  the  sun,  one  after  the  other,  as  the  rocket 
vehicle  spun  in  its  near-vertical  trajectory.  Heavy  lines  indicate 
the  values  of  solar  flux  believed  with  confidence  to  be  important  in 
determining  density  values;  thin  lines  indicate  flux  determinations 
in  which  there  is  lesser  confidence  in  their  importance.  The  diminish- 
ment  of  confidence  is  based  on  several  criteria: 

1.  Slope  of  the  flux  value’s  logarithm  versus  integrated 

2 

number  density  (absorbing  particles,  per  cm  ,  of 
gas  between  the  sun  and  detector)  changes  abruptly, 
indicating  a  sharp  unexplainable  change  of  mean 
absorption  coefficient. 

2.  Flux  values  are  physically  inconsistent. 

3.  Density  determinations  from  the  several  flux 
bands  are  mutually  inconsistent. 

Information  content  of  data  plotted  in  Figs.  1,  2,  and  3  contains  only 
frequencies  below  0.125  Hz. 

Examples  of  how  these  criteria  are  applied  are  revealed  by  the 
curves  in  Figs.  1,  2,  and  3.  A  change  of  slope  in  the  semilog  plot 
of  flux  measured  by  Channel  15  at  about  104  km  is  noted  in  records 
of  experiments  NC7.178  and  NC7.181,  and  at  about  93  km  in  the  record 
of  NC7.180.  In  determining  density  from  the  H  Ly-o  flux  values  at 
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ROCKET  NC  7. 178! INFORMATION  CONTENT  0  I25H 
DETECTOR  B 


FIGURE  2.  Solar  Flux  Versus  Altitude  -  NC7.180  Rocket  Vehicle  Experiment 
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FIGURE3.  Solar  Flux  Versus  Altitude  -  NC7.181  Rocket  Vehicle  Experiment 


higher  altitudes,  an  apparent  bulge  of  02  density  (as  different  from 

the  density  step)  that  33  inconsistent  with  the  determination  of 

"total”  particle  density  (i.e.,  sum  of  number  densities  of  0 and 

N2  and  half  the  number  density  of  0)  occurs.  This  apparent  bulge  is 

reported  by  others  (Lockey,  Horton  &  Rofe,  1969).  In  the  case  of  the 

present  experiments  it  is  believed  to  be  the  effect  of  absorption 

contamination  of  the  H  Ly-e*  line.  This  line  penetrates  deeply  into 

the  atmosphere  because  of  its  small  absorption  cross  section  (lx  10" 
2 

cm  ).  Other  radiation  within  the  sensitivity  band  of  the  H  Ly-cy  de¬ 
tector,  e.g.,  the  solar  Si  III  line  at  1206. 5A,  may  be  the  source  of 
the  concamination.  The  latter  line,  possessing  perhaps  2  percent  of 
the  energy  (Hinteregger,  Hall  6-  Schmidtke,  1964)  of  the  H  Ly-c*  line, 
has  an  absorption  cross  section  2000  times  larger  (Watanabe,  1958). 
Although  solar  energy  of  the  1206. 5A  line  disappears  quickly  as  radi¬ 
ation  penetrates  the  absorbing  atmosphere,  its  effects  perhaps  con¬ 
taminate  measurements  of  the  H  Ly-a  line  by  the  H  Ly-(y  detector. 

Another  example  of  explainable  inconsistency  in  flux  appears  in 

the  record  of  SR  continuum  flux  measured  by  NC7.181  Channel  18 

Detector  B,  shown  in  Fig.  3.  In  this  record,  the  flux  diminishes  by 

-2  -1 

absorption  from  its  out-of-the-atmosphere  value  of  1  erg  cm  sec 

-4  -2  -1 

at  about  150  km,  to  about  10  erg  cm  sec  at  about  106  km,  and 

_2 

then  increases  as  it  penetrates  to  a  local  maximum  of  about  2(10  ) 

-2  -1 

erg  cm  sec  at  approximately  90  km.  This  anomalous  "negative 
absorption"  is  due  to  airglow-like  re-emission  of  SR  continuum  radia¬ 
tion  by  a  third-body  recombination:  atomic  oxygen  created  by  the 
H  Ly-e*  produced  dissociation  of  0^  molecules  recombines  to  form 
molecular  oxygen.  Similar  re-emissions  of  SR  continuum  radiation 
at  altitudes  below  about  100  km  are  apparent  in  the  records  of 
NC7.178  and  NC7.180. 

B.  SOLAR  FLUX  OUTSIDE  THE  ABSORBING  ATMOSPHERE 

Solar  flux  intensities  measured  outside  the  absorbing  atmos¬ 
phere  are  given  in  Table  2.  Measurements  similar  to  NC7.178,  180, 

&  181  data  in  the  1-20^  wave  band  are  given  in  Fig.  4.  This  figure 
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TABLE  2.  FLUX  INTENSITY  OUTSIDE  ATMOSPHERE  ON  21  FEBRUARY  1966 
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presents  solar  flux  measured  by  an  NRL  X-ray  monitoring  satellite 
(1965-93A,  Explorer  30,  NRL  SOLRAD  8)  on  the  day  the  satellite  was 
launched,  which  coincided  with  the  launch  date  of  NC7.178,  180,  &  181. 

Explorer  30  provided  numerous  measurements  of  solar  X-ray  emis¬ 
sions  during  1966  (Kreplin,  1966)  including  those  in  the  1-8,  8-16, 
and  44-60A  wave  bands.  Fig.  4  shows  a  plot  of  X-ray  flux  versus  time 
(UT)  (McClinton  &  Kreplin,  1967)  for  these  wave  bands  on  which  the 
1-20A  measurement  ranges  of  NC7.178,  180  &  181  are  superimposed  for 
time-of -measurement  comparison.  The  ratio  of  1-20A  to  O.sK  flux  is 
0.05  at  NC7.181  launch  time  and  0.035  at  NC7.178  and  NC7.180  launch 
times.  Flux  was  essentially  constant  between  1300  and  2000  GMT.  This 
period  includes  the  launch  times  of  NC7.178,  180,  and  181.  Other  data 
establishing  solar  emission  constancy  during  the  experimental  period 
are  given  in  Appendix  A  and  by  Grobecker  (1967d). 

NC7.178,  180,  and  181  data  provide  measurements  of  1-20^  flux 

that  are  comparable  with  those  of  Explorer  30  when  computed  for  a 

6  o 

sun  temperature  of  2  x  10  K  in  this  band.  Computation  of  flux  at 
a  higher  sun  tempe  -\ature  (4  x  106  °K)  makes  for  better  consistency 
of  density  data  at  lower  altitudes.  Consequently,  at  the  suggestion 
of  Chubb  (1970),  flux  in  the  1-20)5  band  is  computed  as  though  it 
originated  from  two  sources  using  a  certain  ratio  of  out-of-the- 
atmosphere  emission  intensities.  This  ratio  is  the  relationship  be¬ 
tween  intensities  measured  at  the  same  time  by  Explorer  30  and  given 
in  Fig.  4,  and  intensities  with  emission  temperatures  of  2  x  10b  °K 
in  the  8-30^5  band  and  of  4  x  10 6  °K  in  the  1-8A  band. 

For  another  comparison  of  X-ray  data,  a  plot  prepared  from  data 
by  Kreplin  &  Gregory  (1963)  of  X-ray  solar  flux  measurements  by 
other  experimenters  is  shown  in  Fig.  5.  Observations  of  2-8^  radi¬ 
ation  intensity  vary  over  5  orders  of  magnitude  in  Fig.  5,  and  those 
of  0.5-3  and  1-20A  vary  over  nearly  3  orders  of  magnitude.  Computa¬ 
tion  of  these  radiations,  as  being  from  a  gray  body  source,  are  sen¬ 
sitive  to  the  assumption  of  effective  temperature  of  one  or  more 
sources  in  the  sun. 
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C.  COMPUTATION  OF  SOLAR  FLUX 

Reduction  of  measurement  data--an  involved  process  accounting 
for  laboratory  calibrations  of  detectors  and  amplifiers  and  the  ef¬ 
fect  of  solar  aspect  angle  with  respect  to  detector  axis — is  de¬ 
scribed  in  detail  elsewhere  (Grobecker,  1967b,  1967c,  and  1967e). 
Results  of  the  measurements,  expressed  as  detector  current  corrected 
to  calibration  conditions  for  sensitivity  to  normally  directed  flux, 
are  summarized  in  the  plots  of  Figs.  6,  7,  and  8.  Time  is  plotted 
against  detector  current  (at  an  8-cps  rate)  which  is  corrected  for 
detector  aperture  function  and  solar  aspect  angle  as  measured  in 
each  of  two  detectors  for  each  of  the  four  range  bands  measured. 

Each  plot  point  in  Figs.  6  to  8  represents  the  mean  of  20  or  more 
points  selected  from  continuously  recorded  data,  digitized  at  about 
400  cps.  Since  the  44-60/J  flux  determinations  proved  disappointingly 
sensitive  to  errors  in  determination  of  1-20A  flux,  the  results  from 
the  44-60A  detectors  (as  prepared  for  these  experiments)  are  believed 
not  useful  for  determining  44-60^  flux  and  are  not  discussed  in  the 
remainder  of  this  report. 

Corrected  detector  current  data  of  Figs.  6,  7,  and  8  are  used 
in  evaluating  solar  flux,  measured  by  the  three  experiments,  as  a 
function  of  time.  In  the  case  of  the  1-20^  wave  band,  the  evaluation 
of  solar  flux  involves  the  application  of  Eq.  1: 


where 


I 

T, 


=  ion  chamber  current  due  to  radiation  in  1-30^  band. 
=  temperature  of  solar  source  of  1-8^  radiation  = 

4  x  106  °K 
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FIGURE  6.  Detector  Current  Versus  Time  -  NC7 . 178 
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FIGURE  7.  Detector  Current  Versus  Time  -  NC7.180 


CORRECTED  CURRENT  00* 


lYMAN-o  DETECTOR 


T,n  =  temperature  of  solar  source  of  8-30A  radiation  = 

2  x  106  °K 

6(X)  =  efficiency  of  "l-20i  detector”  of  normalized  spectral 

shape  given  in  Fig.  9  and  value  at  8A  given  in  Table  3. 

q  =  charge  on  the  electron 

c  =  1.43879  cm  °K 

10 

w  =  5.4(10)  ion  pairs  per  erg  based  on  ionization 
potential  11.47V  for  filling  gas  CCl^ 
a  =  area  of  detector  aperture  as  given  in  Table  3. 

Gg  =  gas  gain,  as  given  in  Table  3. 

G  =  gain  reduction  due  to  aperture  restriction,  as  given 

in  Table  3. 

o(X)  =  0.22an  (\)  +  0.78aw  (x)  »  determined  from  values  given 
2  w2 
in  Fig.  10. 

N(z^)  =  T  FCz^JnCz^dz^  -  overhead  air  mass,  above  altitude  z^; 
Zi 

based  on  U.S.  Standard  Atmosphere,  1962. 

FCz^  =  Optical  depth  factor  given  in  Table  4. 

K  =  E(l-8A)(z^)/E(8-30i)(z^),  ratio  of  flux  in  1-8A  and 

8-30^  bands  as  determined  from  measurements  of  Explorer 
30  satellite,  given  in  Fig.  4. 

Evaluation  of  solar  flux  in  the  SR  continuum  involves  the  appli¬ 
cation  of  a  relation  similar  to  that  of  Eq.  1,  except  the  solar 
source  is  taken  to  be  single  and  at  a  temperature  of  4800  °K.  The 
spectra  of  detector  efficiency  for  the  detectors  of  the  SR  continuum 
are  given  in  Fig.  11. 

H  Ly-o  flux  is  evaluated  in  terms  of  ion  chamber  current  by: 


I  =  qwa€E 


(2) 


where  I  =  ion  chamber  current 

€  =  photon  detection  efficiency  at  H  Ly-a 
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TABLE  3.  EXPERIMENT  INSTALLATION  PARAMETERS 
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TABLE  4.  OPTICAL  DEPTH  FACTORS 
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E  =  energy  flux  in  Lyman-a  line 

w  =  ion  pairs  per  erg,  assumed  to  be  one  per  photon 
a  =  window  area,  accounting  for  aperture  restrictions 
q  =  electron  charge. 

The  spectrum  of  detection  efficiency  of  the  H  Ly-a  detectors  is  given 
in  Fig.  12. 

These  computations  yield  data  of  point  frequency  about  8  cps, 
which  ait;  ncisy  and  require  smoothing.  The  data  were  smoothed  by  a 
sequence  of  filtering  operations:  frequencies  of  the  measurement 
spectrum  expected  to  be  artifacts  of  the  instrumentation  are  elimi¬ 
nated,  and  frequencies  expected  to  be  characteristic  of  the  geophysi¬ 
cal  variation  are  left  intact. 

The  progressive  elimination  of  unwanted  frequencies  from  the 
data  is  illustrated  in  the  frequency  chart  of  Fig.  13.  The  maximum 
frequency  content  of  the  original  measurement  data  is  about  200  cps, 
corresponding  to  the  approximate  5  ms  rise  time  of  the  detector.  On 
this  basis,  the  continuous  record  of  telemetered  data  is  digitized 
at  400  cps  to  ensure  that  all  frequencies  contained  in  the  instrument 
signal  are  contained  in  the  digitized  data. 

Spectral  analysis  of  the  digitized  data,  of  which  Fig.  14  is  an 
out-of-the-atmosphere  example,  shows  that  frequency  power  greater 
than  50  cps  in  the  measurement  data  is  negligible.  Spectral  analyses 
of  data  sections  from  each  rocket  digitized  at  800  cps  showed  no 
additional  lines.  It  is  believed  that  no  "aliasing"  frequencies 
exist  in  the  spectrum  of  data  digitized  at  400  cps. 

Frequencies  near  50  cps,  shown  in  Fig.  14  and  schematically  in 
Fig.  13,  are  harmonics  ascribed  to  irregularities  in  the  ion-chamber 
aperture  function.  This  function  describes  the  reduced  sensitivity 
of  the  detector  to  radiation  incident  at  angles  other  than  the  normal 
incidence  ones  of  laboratory  calibrations,  and  its  variations  depend 
on  detector  geometry  and  composition,  which  must  be  accounted  for 
(Grobecker,  1967e). 
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ATMOSPHERE  CHARACTERISTICS  f  -g-cpt 

FREQUENCY  OF  ROCKET  SPIN  8  cpt 
SPECTRAL  ANALYSIS  SHOWS  f  50  cp, 
FROM  APERTURE  FUNCTION 
DIGITIZED  AT  f  *  400  cp* 

RISE  TIME  OF  DETECTOR  5  rm 


INITIAL  FREQUENCIES  HAVE  HIGH  FREQUENCY  CUTOFF 
I*  200  BY  INSTRUMENTATION 


APERTURE  FUNCTION  CORRECTION  PRODUCED  DATA 
CONTAINING  ONLY  f  •:  8  cps  IN  DATA  AT  f  «  8  cpt 


"HAMMING  8"  AND  "DECIMATION  8Y  4"  PRODUCES  DATA 
CONTAINING  ONLY  f  <  1  cpt  IN  DATA  AT  f  *  2  epi 


"HAMMING  2"  AND  "DECIMATION  8Y  2"  PRODUCES  DATA 
CONTAINING  ONLY  f  0.5  cpi  IN  DATA  AT  f  *  1  cpt 


(^aI^^aI 


"HAMMING  2"  AND  "DECIMATION  BY  2"  PRODUCES  DATA 
CONTAINING  ONLY  f  -  0.25  cpi  IN  DATA  AT  f  =  0.5  cpt 


FIGURE  13.  Smoothing  of  Solar  Flux  E(t) 
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NORMALIZED  SPECTRAL  DENSITY  ESTIMATES 


1.25 


FREQUENCY,  cpi 


FIGURE  14.  Spectral  Analyses  of  Digitized  Data  of  Rocket  NC7.178 


The  frequency  of  the  rocket  spin  at  about  8  cps  and  its  second 
harmonic  appear  as  a  strong  modulation  of  the  measurement  data.  A 
slow  modulation,  due  to  precession  of  the  rocket  and  causing  yaw 
angle  with  respect  to  the  trajectory,  is  superimposed  on  the  8  cps. 

This  precession  has  a  period  of  more  than  20  sec.  In  addition  to 
the  precession,  nutation  of  the  rocket  (giving  small  yawing  motion) 
is  characterized  by  a  frequency  of  to  cps.  Considering  tne  speed 
with  which  the  rocket  moves  through  the  atmosphere,  and  that  the  ex¬ 
pected  atmospheric  changes  are  to  be  sensed  in  vertical  intervals 
approximated  by  atmospheric  scale  heights,  the  investigated  atmos¬ 
pheric  characteristics  are  expected  to  involve  frequencies  smaller 
than  about  cps. 

To  selectively  filter  out  the  several  unwanted  frequencies  that 
are  artifacts  of  instrumentation  (schematically  shown  in  Fig.  13),  fil¬ 
ters  described  by  Blackman  and  Tukey  (1958)  were  employed  in  a  succession 
of  computing  operations.  These  are  the  "hamming"  filters  D^g  (where 
m  =  8)  and  D^2  (where  m  =  2) .  The  D^Q  filter  converted  320  points  of 
data,  with  content  in  the  0-  to  8 -cps  range,  to  32C  points  of  data 
with  information  content  of  less  than  1  cps  frequency.  Filter  D32 
converted  80  points  of  data  with  content  from  0-1  cps  to  80  points 
with  information  content  of  frequency  less  than  0.5  cps.  Following 
each  hamming,  a  "decimation"  by  4  to  2  cut  down  the  number  of  data 
points.  If  done  properly,  decimation  does  not  introduce  aliased 
frequencies  into  the  information  content  of  the  reduced  band. 


Frequency  transmission  characteristics  of  filters  D^g  and  D^, 
as  applied  to  8  cps  data  and  1  cps  data,  respectively,  are  shown  in 
Fig.  15.  As  evidenced  by  the  Fourier  transforms,  the  content  re¬ 
maining  after  hamming  is  trivial  at  frequencies  higher  than  the  cut¬ 
off  frequency. 

An  N  point  hamming  filter  is  defined  as  follows:  if 


Xj  transforms  by  hamming  to  x^  (i  =  l,...,m), 
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FIGURE  15.  Fourier  Transforms  of  Filters  and  D, 


then  D{  =  0.54  +  0.46  x  cos  £1 


( t  ~  -N ,  N )  and  W 


-£■ 

t"-N 


accordingly  x?  =  x^  (j  =  1,...,N)  and  is  defined  by: 


D  /W  (j  =  N  +  1, . ,.,m  -  N) 


Xj  =  Xj  ^  =  m  "  N  +  *  (3) 

A  decimation  by  N  takes  every  Nth  point: 

Xj  transforms  by  decimation  to  x*  ^  =  1,...,^  ; 

and  x*  =  xx  +  .  1)N  (‘  =  1**-*>S)  * 

A  hamming-N  operation  changes  the  information  content  by  and 
a  decimation  by  N  changes  the  data  frequency  by  Table  5  describes 
the  effects  of  the  successive  filtering  process  applied  to  the  orig¬ 
inal  8  cps  data. 

In  summary,  the  filtering  procedure  applied  to  the  experimental 
data  was  sequenced  as  follows: 

1.  Hamming-8  filter  reduced  the  data  frequency  content 
from  more  than  1  cps  to  less  than  1  cps 

2.  Decimation  by  4  reduced  the  number  of  data  points 
to  be  handled  from  8  to  2  cps 

3.  Hamming-2  filter  reduced  the  data  content  to  fre¬ 
quencies  less  than  0.5  cps 

4.  Decimation  by  2  reduced  the  2  cps  data  to  less 
than  1  cps. 
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TABLE  5.  EFFECTS  OF  SUCCESSIVE  FILTERING  AND  DECIMATION 
OPERATIONS  IN  SMOOTHING  FLUX  DATA  [E(t)] 


Filter  Operation 

Da  ta 

Frequency 

(cps) 

Maximum  (cps) 
Information 
Content 

Hamming-8 

H8 

8 

1 

Decimation  by 

2 

D4 

2 

1 

Hamming-2 

H2 

2 

0  5 

Decimation  by 

2 

D2 

1 

0.5 

Hamming-2 

h2 

1 

0.25 

Decimation  by 

2 

D2 

0.5 

0.25 

Hamming-2 

h2 

0.5 

0.125 

Decimation  by 

2 

D2 

0.25 

0.125 

Hamming -2 

H2 

0.25 

0.062 

Final  data,  whose  frequency  content  was  less  than  0.125  cps  in 
a  data  frequency  of  0.5  cps,  were  produced  by  second  and  third  appli¬ 
cations  of  steps  3  and  4. 

An  example  of  the  data  without  filtering  or  decimations,  digit¬ 
ized  at  8  cps  and  describing  information  of  frequencies  up  to  4  cps, 
is  given  in  Fig.  16.  These  data,  after  hamming  and  decimation,  are 
shown  in  Fig.  17  where  information  content  is  less  than  0.125  cps  at 
a  0.5  cps  data  rate.  Further  iterations  of  hamming  and  decimation, 
which  go  too  far  by  reducing  information  content  to  frequencies  less 
than  0,062  in  data  ooints  of  0.25  cps,  yielded  results  shown  in  Fig. 
18. 

After  reviewing  the  results  of  successive  filtering  applications 
to  all  data,  the  flux  data  of  all  channels  were  reduced  to  data 
points  at  a  frequency  of  0.5  cps,  retaining  information  smoothed  to 
contain  only  frequencies  less  than  0.125  cps. 
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ME  ( *ec ) 


ROCKET  NC7  . 178 


FLUX  («rgs  cm’^wc  S' 


FIGURE  16.  Flux  Versus  Time  at  Data  Frequency  8  cps  and  Information 
Content  1  cps  -  1-20A  Band,  NC7.178  Detector  li, 
Channel  14 
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FLUX  (erg»  cm“^»ec  *) 


FIGURE  17.  Flux  Versus  Time  at  Data  Frequency  0.5  cps  and  Information 
Content  0.1 25  cps  -  1**20A  Band,  NC7.178  Detector  B# 
Channel  14 
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FIGURE  18.  Flux  Versus  Time  at  Data  Frequency  0.25  cps  and  Information 
Content  0.062  cps  -  1-20A  Band,  NC7.178  Detector  B, 
Channel  14 
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D.  FLUX  VERSUS  ALTITUDE 

A  next  step  in  the  determination  is  to  convert  flux,  as  a  func¬ 
tion  of  time  of  rocket  flight,  to  flux  as  a  function  of  altitude, 
using  the  altitude  versus  time  characteristics  given  in  Fig.  19. 

The  altitude  history  of  flux  measured  by  the  three  rocket  ex¬ 
periments  for  the  three  wave  bands  of  interest  are  shown  in  Fig.  1. 
These  flux-altitude  plots  are  derived  from  the  data  of  flux  versus 
time,  smoothed  to  contain  only  information  of  frequency  smaller  than 
0.125  cps. 

Note  that  in  interpreting  density  computations,  in  terms  of 
flux  and  flux  gradient,  the  vertical  velocity  of  the  rocket: 

^  (t)  s  1.3  km  sec”1  (5) 

corresponds  to  a  speed  of  Mach  4;  accordingly,  flux  altitude  data 
are  contained  within  a  band  of  maximum  wave  number  k  ^  0.1  km-1  . 

The  atmospheric  scale  height  implies  atmospheric  variations  are  ex¬ 
pected  within  a  band  of  maximum  wave  number  k  &.  1/6  km-1  .  Rocket 
nutation  contains  wave  numbers  as  small  as  k  *  1/4  km-1  .  All  other 
undesirable  wave  numbers  (k  >  1/4  km”1)  are  those  associated  with 
rocket-  and  instrumentation-induced  characteristics .  These  spectral 
differences  provide  a  basis  for  distinguishing  instrumentation  arti¬ 
facts  from  geophysical  effects. 

Some  editing  of  flux  versus  altitude  curves  appears  desirable  be¬ 
fore  beginning  density  computations.  In  particular,  all  flux  and 
flux  values  smaller  than  10”1(1  are  ignored.  By  accounting  for  all 
factors  previously  described,  including  computation  results  of  vari¬ 
ation  coefficients  described  in  the  following  paragraphs,  credible 
ranges  of  measurement  for  each  detector  can  be  derived  (Table  6). 

Data  plotted  in  heavy  lines  in  Figs.  1  to  3  inclusive,  within  the 
credible  range  limits  described  in  Table  6,  are  deemed  reliable  for 
computing  number  density. 


35 


TABLE  6.  CREDIBLE  RANGES  OF  ALTITUDE  AND  TIME 


Rocket 

Channel 

Detector 

— 

Plot 

Limit 

Low  Alt. 
Limit 
(Z) (km) 

Upper  Alt. 
Limit 
(Z) 

Low  Time 
Limit 
(t) 

Upper  Time 
Limit 
(t)(sec ) 

NC7.178 

189 

A 

103 

123 

68.2 

84.2 

Kfl 

B 

107 

123 

71.3 

84.5 

1 

A 

70 

97 

45.1 

63.7 

97* 

104* 

63.7* 

69.0* 

B 

70 

100 

43.1 

66.0 

100* 

110* 

66.0* 

73.6* 

18 

.  A 

115 

129 

77.6 

89.5 

B 

122 

141 

83.4 

100.7 

NC7.180 

14 

A 

98 

121 

66.7 

85.8 

B 

98 

121 

66.7 

85.8 

15 

A 

70 

93 

46.5 

62.9 

93* 

105* 

63.0* 

72.3* 

B 

70 

95 

46.5 

64.4 

95* 

105* 

64 . 4* 

72.3* 

18 

A 

108 

140 

74.7 

104.4 

B 

108 

147 

74.7 

112.1 

NC7.181 

14 

A 

98 

125 

67.5 

90.6 

B 

98 

125 

67.5 

90.6 

15 

A 

85 

96 

57.7 

66.0 

96* 

110* 

66.0* 

77 . 3* 

B 

86 

99 

58.4 

68.3 

99* 

110* 

68.3* 

77.3* 

18 

A 

125 

141 

91.0 

107.0 

B 

117 

.i.44 

83.1 

110.4 

*Reduced  confidence  -  see  Section  III. A.  D  B. 
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E. 


ERROR  ESTIMATES 


An  analysis  of  estimated  error  contributions  (fractional  uncer¬ 
tainties)  to  the  measurement  of  flux  by  each  detector  of  each  channel 
is  summarized  in  Table  7. 

In  the  analysis,  data  variation  coefficients  for  each  detector 
used  or  for  available  representative  samples  of  like  kind  are  re¬ 
corded,  and  they  comprise  several  types  including  systematic  varia¬ 
tions  in  laboratory  techniques  for  determining  detector  sensitivity 
and  variations  in: 

•  Sensitivities  of  individual  detectors 

•  Gas  gain  in  the  time  between  laboratory  calibration 
and  time  of  flight 

•  Detector  aperture  area  and  shape  of  aperture  function 

•  Gain  in  dc  amplifiers 

•  Telemetry  timing  and  other  variations  introduced  by 
telemetry  and  digitizing  processing. 

Based  on  the  premise  that  each  of  these  contributions  is  inde¬ 
pendent  of  the  others,  the  total  variation  coefficient  of  normalized 
ion  current  is  computed  from  the  rms  of  all  contributions,  for  com¬ 
parison  with  the  measured  variation  coefficient  of  th^  current  values 
determined  from  telemetered  data.  Results  of  the  computed  and  meas¬ 
ured  variation  coefficients  of  the  normalized  detector  current  are 
comparable  in  the  case  of  each  detector. 

Variation  coefficients  of  flux  values  over  the  credible  range 
as  determined  by  each  single  detector  are  also  tabulated  in  Table  7. 
Differences  between  these  variation  coefficients  and  those  of  detec¬ 
tor  current  reflect  distortion  introduced  by  smoothing  and  filtering 
processes  in  time. 

The  variation  coefficient  of  measured  flux  data  was  determined 
in  the  following  way. 

The  coefficient  of  variation,  or  fractional  variation  (i.e., 
the  ratio  of  standard  deviation  to  the  mean)  in  the  filtered  flux 
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TABLE  7.  ESTIMATED  FRACTIONAL  UNCERTAINTIES  OF  FLUX 
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1.0  0.06  (I.r  0  0.02  0.050  0  0.007  1.0  0.030  0.20 


data  is  estimated  as  a  function  of  time  for  rocket  experiments 
NC7.178  in  Fig.  20,  NC7.180  in  Fig.  21,  and  NC7.181  in  Fig.  22.  The 
error,  as  shown,  is  computed  for  8-sec  intervals  within  an  altitude 
range  near  the  credible  ranges  of  altitude  shown  in  Table  6. 

Computation  of  fractional  variation  is  based  on  a  calculation, 
for  each  time  increment  T  ,  of  the  difference  between  unfiltered 
flux  values  (XUF^)  and  filtered  flux  values  after  hamming  and  deci¬ 
mations  (XF^),  wherein  t  takes  values  from  1  to  NP  where  NP  is  the 
number  of  points  in  the  rocket  data;  N  is  determined  for  each  T  so 

i 

that: 


T  „  £  T  +  8  <;  T 
t  +  N  t  t 


+  N  +  1 


(6) 


and  the  error  is  calculated  as: 


(XUF..  -  XFj  )2/ (N  +  1) 


XFj/(N  +  1) 


(7) 


In  addition  to  the  fractional  variation  computed  by  the  rela¬ 
tion  (7),  a  single  coefficient  of  variation  is  calculated  for  the 
entire  credible  range. 


40 


FRACTIONAL  ERROR  OF  FLUX _  _ FRACTIONAL  ERROR  OF  FLUX  FRACTIONAL  ERROR  OF  FLUX 
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FIGURE  20.  NC7.178  Flux  Error  Versus  Time 


FRACTIONAL  ERROR  OF  FLUX  _  _  _ FRACTIONAL  ZRROR  OF  FLUX  FRACTIONAL  ERROR  OF  FLUX 


FIGURE  21 .  NC7. 180  Flux  Error  vs.  Time 
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FIGURE  22.  NC7.181  Flux  Error  vs.  Time 
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III.  "TOTAL"  PARTICLE  AND  02  DENSITY  DISTRIBUTION 
A.  ANALYSIS  OF  DATA 

Reduction  of  the  solar  flux  absorption  data  is  described  in  de¬ 
tail  by  Grobecker  (1967b,  1967c,  and  1967e).  Briefly,  the  telemetry 
data  of  ion  chamber  signals  continuously  recorded  at  the  ground  sta¬ 
tion  during  the  rocket  flight  were  reduced  in  the  following  steps: 

1.  Digitized  at  400  samples  per  second  in  terms  of  volt¬ 
age  of  the  rocketborne  dc  amplifiers 

2.  Converted  by  amplifier  calibrations  to  ion  chamber 
detector  output  current 

3.  Normalized  to  equivalent  flux  energy  in  the  wave¬ 
length  band  of  detector  sensitivities,  accounting 
for  laboratory  calibrations  of  detector  current 
resulting  from  normally  incident  flux,  variation 
of  detector  sensitivity  when  the  flux  is  not 
normally  incident,  and  the  time  history  of 
individual  detector  changes . 

Determination  of  equivalent  flux  energy  also  required  the  calculation 
of  the  effect  of  solar  radiation,  progressively  hardening  (mean  wave¬ 
length  within  the  band  of  detector  sensitivity  becomes  shorter)  as  it 
penetrates  deeper  into  the  atmosphere,  on  the  apparent  spectral  effi¬ 
ciency  determined  by  laboratory  calibration  for  each  individual  detector. 

From  data  on  flux  E(z)  as  a  function  of  altitude  (Section  II.D.), 
density  of  the  "total"  particle  and  0^  content  of  the  atmosphere  may 
be  computed.  In  this  case,  "total"  particle  content  is  taken  to  be: 

n(tot)  =  n(N2)  +  n(02)  +  Jg  n(0)  ,  (8) 
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since  the  X-radiation  does  not  distinguish  between  the  atomic  and 
molecular  form  of  oxygen.  Density  is  computed  by  relation: 


n(z) 


dE/dz 

ffC'z)  F(z)  E(z) 


(9) 


where  a(z)  is  the  mean  absorption  cross  section  as  a  function  of  al¬ 
titude,  F(z)  is  the  optical  depth  factor,  and  E  is  the  solar  flux 
intensity. 

The  evaluation  of  dE/dz(z)  from  the  digitized  values  of  flux 
E(z)  again  required  careful  band-limited  filtering  accomplished  in 
the  wave  number  domain  by  means  of  fast  Fourier  transforms  (Cooley  & 
Tukey,  1965).  The  band  of  wave  numbers  retained  in  the  computation 
of  dE/dz(z)  included  only  those  between  0  and  1/7  cycles  per  kilo¬ 
meter.  In  computing  density  in  terms  of  flux  and  flux  gradient, 
advantage  is  taken  of  the  fact  that  major  atmospheric  characteristics 
are  expected  to  be  observed  within  an  atmospheric  scale  height  of 
7  km  to  10  km  and,  hence,  are  contained  within  a  band  of  wave  numbers 
lers  than  k  =  1/7  cp  km.  Rocket  nutation  is  characterized  by  wave 
numbers  larger  than  k  a  1/4  cp  km.  All  other  undesirable  wave  num¬ 
bers,  k  >  1/4  cp  km,  are  taken  to  be  those  associated  with  rocket- 
and  instrumentation- induced  characteristics,  and  not  geophysical 
ones.  Data  filtering  in  terms  of  these  spectral  differences  provides 
a  basis  for  distinguishing  instrumentation  artifacts  from  major  geo¬ 
physical  effects,  but  it  obscures  observation  of  the  fine  structure 
of  the  atmosphere. 

An  essential  in  computing  density  distribution  by  Eq.  9  is 
the  appropriate  evaluation  of  the  optical  depth  factor  F(z).  This 
factor  approximates  the  secant  of  the  solar  zenith  angle  at  small 
values  of  the  angle,  but  also  accounts  for  the  apparent  thickness  of 
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the  spherically  shelled  atmosphere  at  large  zenith  angles.  Optical 
depth  factors  were  calculated  from  tables  and  with  methods  descrj  bed 
by  Swider  (1964),  accounting  for  local  solar  zenith  angles  at  times 
of  the  rocket  flight,  and  an  assumed  model  of  atmospheric  scale 
height  H(z). 

The  mean  absorption  cross  section,  except  for  that  of  the  H  Ly-<y 

-20  2 

line  which  is  taken  at  a  constant  value  of  1  x  10  cm  over  the 
width  of  the  line  (Watanabe,  1958),  is  calculated  by  the  relation: 


-o(\)  /  F(z)  n(z)  dz  ;  -'j(X)  f  F(z)  n(z)  dz 

k  f  P(\,Ta)  o(\)e  zi+l  dx  +  f  P(X,TJ0)  i(X)e  zi  +  l  dx 

'(Zj)  =  -hi - - - * - 

/  F(z)  n(z)  dz  '  ..(X)  /  F( z )  n(z)  az 

k  /  P(X,T8)e  zi  +  l  dX  ♦  J  P(X,TJ0)  o(X)e  Zi  +  1  dx 

IA  U* 

(10) 

where  k  =  E(1-8A)/E(8-30A) ,  the  ratio  of  out-of-the- 

atmosphere  fluxes  of  the  1-8A  and  8-30A  wave 
bands,  as  determined  for  the  same  hour  by 
Explorer  30  satellite. 

=  0.05  for  NC7.181  and  0.0035  for  NC7.178  & 

NC7.180. 

P(X,T8)  =  \"5  [exp(c/T8\)-l]'1 

P(X,T30)  =  X"5  Ccxp(c/T30X)-l]"1 

Tg  =  2(10)6  °K 

T30  =  4(10)6  °K 

c  =  1.43879  °K 
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T  =  solar  temperature 

c^X)  =  wavelength-dependent  absorption  cross  section. 

CO 

and  J  n(z)  dz  is  derived  from  values  of  density  in  previous 
zi+l 

(higher  altitude)  computational  steps.  Relation  10  expresses  recogni¬ 
tion  that  the  mean  wavelength  within  the  band  of  detector  sensitivity 
becomes  shorter  as  the  solar  radiation  penetrates  deeper  into  the 
atmosphere  (hardening). 

For  1-20 A  determinations,  solar  temperature  T  is  specified  for 

#  •  fi 

two  separate  sources:  1-8A  at  Tg  =  8  x  10  and  8-30A  at  =  2  x  10 
°K. 

For  the  1-20A  radiation,  the  effective  absorption  cross  section 
is  due  to  the  effect  of  both  02  and  N2  molecules  and  to  the  effect  of 
0  atoms.  The  cross  section  of  0  atoms  is  half  that  of  02  molecules. 
Therefore, 


a(\)  =  0.22  (X)  +  0.78  ctn  (X) 


(ID 


The  wavelength-dependent  cross  section  c(X),  used  in  determining 
mean  cross  section  cr(z)  by  means  of  Eq.  10,  is  given  in  Table  8 
utilizing  data  of  Hinteregger  et  al.  (1964). 

For  the  SR  continuum,  values  of  wavelength-dependent  absorption 
cross  section  a(X)  given  in  Table  9  (Hinteregger  et  al. ,  1964)  are 
used  to  compute  the  mean  absorption  cross  section  a(X)  over  the  band 
of  detector  sensitivity  for  radiation  in  the  band  of  a  solar  source 
at  temperature  4800  °K.  This  computation  uses  an  equation  analogous 
to  Eq.  10  with  a  single  term  on  the  right-hand  side  numerator  and 
denominator. 
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TABLE  8.  ABSORPTION  CROSS  SECTION  a(\)  IN  1-2oA  RANGE 


\(A) 

an  (cm2) 
u2 

aM  (cm2) 

n2 

ct(\) 

1-3 

0.002  x  10" 18 

0.0015  x  10" 18 

0.0015  x  10"18 

3-5 

0.012  x  10" 18 

0.0065  x  10" 18 

0.0075  x  10" 18 

5-10 

0.075  x  10" 18 

0.045  x  10“18 

0.040  x  10"18 

10  -  15 

0.27  x  10"18 

0.15  x  10"18 

0.171  x  10"18 

15  -  22 

0.7  x  10"18 

0,36  x  10-18 

0.42  x  10" 18 

22  -  31 

0.09  x  10"18 

1.0  x  10"18 

0.80  x  10“ 18 

TABLE  9.  ABSORPTION  CROSS  SECTION  a( X)  IN  SR  CONTINUUM 


\ 

(A) 

a(X) 

(cm2 ) 

1325  - 

1375 

5.8 

X 

H 

O 

l 

H 

en 

1375  - 

1475 

14.3 

X 

ID"18 

1475  - 

152  5 

11.3 

X 

10- 18 

As  in  the  case  of  flux  values  depicted  in  Figs.  1,  2,  and  3  ('Sec¬ 
tion  II),  density  determinations  have  altitude  ranges  of  best  and 
reduced  credibility.  Diminishment,  from  the  best  credibility  shown 
by  heavy  lines  in  Figs.  1,  2,  and  3,  is  based  on  one  of  several  con¬ 
siderations: 

1.  Slope  of  the  flux  value’s  logarithm  versus  integrated 

2 

number  density  (absorbing  particles,  per  cm  ,  of 
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gas  between  sun  and  detector)  changes  rapidly, 
indicating  an  unexplainable  change  of  mean 
absorption  coefficient. 

2.  Flux  values  are  physically  inconsistent. 

3.  Density  determinations  from  the  several  flux 
bands  are  mutually  inconsistent. 

Altitude  ranges  for  which  measurements  of  density  are  credible 
(Table  6,  Section  II. D.)  are  based  on  these  considerations.  Those 
which  reflect  reduced  confidence  are  marked  with  an  asterisk  in 
Table  6. 

B.  DENSITY  VERSUS  ALTITUDE 
1.  Absorption  Measurements 

Values  of  density  computed,  with  Eq.  9,  as  a  function  of 
altitude  are  shown  for  NC7.178,  NC7.180,  and  NC7.181  in  Figs.  23, 

24,  and  25,  respectively.  Data  points  for  each  of  two  detectors  (A 
and  B),  for  each  wave  band  measured  by  the  three  experiments,  are 
indicated  separately  by  symbols.  Also  shown  in  the  figures,  by  hori¬ 
zontal  bars,  are  the  limits  of  the  range  of  credible  data.  These 
limits  are  within  the  wider  range  of  data  that  were  analyzed  computa¬ 
tionally  but  included  analyzed  points  in  adjacent  ranges  believed, 
for  reasons  given  in  this  report,  not  representative  of  the  geophysi¬ 
cal  domain.  Representations  of  density,  taken  as  a  mean  of  the  two 
detector  measurements  in  each  band,  are  indicated  by  thin  lines  in 
the  figures. 

Examples  of  the  application  of  the  criteria  for  credible  ranges 
listed  in  Table  6  (Section  II)  are  found  in  the  flux  values  and  semi- 
log  plots  of  flux  versus  integrated  number  density  in  Figs.  1,  2, 
and  3  and  the  density  values  of  Figs.  23,  24,  and  25.  As  mentioned 
earlier  (Section  II. A, ),  change  of  slope  in  the  semilog  plot  of  flux 
measured  by  Channel  15  at  particular  altitudes  was  noted  in  records; 
also,  in  determining  density  from  H  Ly-a  flux  values,  an  apparent 
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FIGURE  24.  Density  versus  Altitude  -  NC7.180 
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DENSITY  (CM-J) 


FIGURE  25.  Density  versus  Altitude  -  NC7.181 


bulge  of  density  occurred  at  high  altitudes.  As  stated  in  Section 
II. A.,  this  bulge  is  thought  to  be  the  effect  of  absorption  contam'n- 
ation  by  other  radiation  (e.g.,  solar  Si  III  line  at  1206. SA)  in  the 
H  Ly-<y  line  which  penetrates  deeply  into  the  atmosphere.  The  1206.5)5 
solar  energy  disappears  quickly  as  radiation  penetrates  the  absorbing 
atmosphere,  but  its  effect  perhaps  contaminates  measurements  of  the 
H  Ly-cr  line  by  the  H  Ly-a  detector  down  to  the  lower  altitudes  indi¬ 
cated  by  asterisks  in  Table  6  (Section  II). 

2.  Ionosonde  Measurements 

Ionosondes  were  operated  at  Fort  Churchil]  and  at  Eglin  on  the 
day  of  the  launch.  Data  were  analyzed  by  the  ESSA  Research  Labora¬ 
tories,  Boulder,  Colorado  (Wright,  1966).  A  plot  of  the  ionosonde 
electron-density  distribution  at  the  time  and  place  of  the  three 
launches  is  shown  in  Fig.  26. 

Note  that  the  gradient  of  electron  density  increases  greatly 
at  altitude  118  km  for  NC7.178,  107  km  for  NC7.180,  and  108  km  for 
NC7.181. 


C.  GEOPHYSICAL  INTERPRETATION 


The  geophysical  interpretation  of  data  from  the  three  rocketborne 
experiments  is  described  here.  Present  data  are  compared  with  data 
on  neutral  particle  densities  determined  by  other  experimenters,  with 
electron  density  profiles  determined  by  ionosonde  in  the  same  alti¬ 
tude  region,  and  with  other  reoresentative  observations. 


1.  "Total11  Particle  Density  Profile 

Profiles  of  "total"  particle  density  for  each  of  the  three 
rocketborne  experiments  are  shown  in  Fig.  27.  They  are  depicted  in 
comparison  with  the  U.S.  Standard  Atmosphere  (1962)  and  COSPAR  Inter¬ 
national  Reference  Atmosphere  CIRA  (1965),  using  models  appropriate 
to  conditions  for  the  launch  day  (10.6  cm  radiation  flux  of  about 

-  4).  Also  included 


-22  -2  -1 
85  x  10  w  m  Hz  and  geomagnetic  index  A 


in  Fig.  27,  for  comparison  purposes,  are  data  reported  from  mass 
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FIGURE  26  lonosonde  Electron  Density  Profiles — Fort  Churchill  and  Eglin 
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FIGURE  27.  Comparison  of  "Total"  Particle  (number)  Density  with  other  Theoretical 
and  Experimental  Determinations 


spectrometer  measurements  made  at  Fort  Churchill  on  12  December  1966 
(Gross,  Offermann  and  von  Zahn,  1967),  and  in  Sardinia  on  11  December 
1965  (Mauersberger,  Muller,  Offermann,  &  von  Zahn,  1965),  and  data  re¬ 
ported  by  von  Zahn  and  Gross  (1968).  The  "total"  particle  density 
[n(Ng)  +  n(O^)  +  %  n(0)]  profiles  are  based  on  the  assumption  that 
the  solar  source  of  X  rays  in  the  1-2oA  band  is  at  a  temperature  of 
2  x  106  °K  and  the  solar  source  of  X  rays  in  the  1-20^  band  is  at  a 
temperature  of  8  x  10^  °K.  Flux  data  for  density  computations  were 
filtered  to  include  only  frequencies  smaller  than  0.125  cps.  Data  of 
greatest  credibility  are  shown  with  a  heavy  line;  those  of  lesser 
credibility  with  a  thin  line. 

As  evidenced  in  Fig.  27,  determinations  of  "total"  particle  den¬ 
sity  for  the  three  rocketbome  experiments  at  about  120  km  are  con¬ 
sistent  with  each  other  and  with  the  CIRA  1965  determination.  Density 
values  at  lower  altitudes  are  higher  for  the  present  experiments  than 
those  attributed  by  the  U.S.  Standard  Atmosphere  (1962).  But  they 
are  consistent  with  the  high  temperature  from  70  to  125  km,  implied 
by  the  smaller  slopes  d  In  n/dz  displayed  with  mutual  consistency  by 
the  profiles  of  both  "total"  and  0 2  density  below  100  km. 

2.  0 ^  Particle  Density  Profiles 

Profiles  of  02  particle  density  for  each  of  the  three  rocketborne 
experiments  are  shown  in  Fig.  28.  However,  02  density  determinations 
derived  by  two  SR  detectors  of  NC7.178  and  by  one  SR  detector  of 
NC7.181  are  omitted  from  Fig.  28.  These  determinations,  shown  in 
Figs.  23  and  25,  appear  dubious  because  of  the  small  signal  strength 
recorded  by  the  detectors  in  question. 

Heavy  lines  in  Fig.  28  denote  experimental  data  of  best  credi¬ 
bility;  thin  lines  denote  data  of  diminished  credibility.  Horizontal 
bars  demark  the  credibility  ranges. 

Since  the  absorption  of  H  Ly-o  radiation  may  be  contaminated  by 
absorption  of  the  Si  III  line  at  1206. sA,  the  densities  are  suspect 
when  determined  from  H  Ly-»  absorption  above  the  104,  105,  and  110  km 
altitudes  (and  possible  above  97,  93,  and  96  km)  for  NC7.178,  NC7.180, 
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FIGURE  28.  Comparison  of  O2  Number  Density  with  Other  Theoretical  and 
Experimental  Determinations 
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and  NC7.181,  respectively.  This  contamination  by  the  Si  III  line  (2 
percent  of  the  out-of-the-atmosphere  energy)  may  account  for  the 
bulge  of  density  appearing  in  the  H  Ly-a  derived  profiles  of  NC7.178 
and  NC7.181. 

Despite  the  questionability  of  the  bulge  and  the  limitation  of 
the  dynamic  range  of  the  duplicated  detection  systems,  which  pre¬ 
vented  each  detector  type  from  covering  the  whole  range  of  the  02 
density  step,  the  H  Ly-a  detectors  indicate  a  steep  gradient  of  0 2 
density  between  108  and  117  km;  this  is  well  within  the  dynamic  range 
of  density  determination  and  is  reasonably  consistent  with  the  SR 
determination  of  density.  In  the  case  of  NC7.181,  the  overlap  be¬ 
tween  the  minimum  limit  of  density  determinable  by  H  Ly-a  absorption 
and  maximum  limit  of  density  determinable  by  SR  absorption  may  be 
attributed  to  a  lag  introduced  by  the  filtering  process.  Although 
the  existence  of  the  step  in  02  density  is  apparent  from  the  data, 
its  altitude  is  uncertain  within  perhaps  4  km. 

Also  shown  in  Fig.  28,  for  comparison  purposes,  are  data  points 
reported  by  other  observers  (Lockey  et  al. ,  1969  and  Quessette,  1970), 

and  standard  atmosphere  models  [CIRA  1965  model  appropriate  to  a 

-2  -1 

10.6-cm  flux  of  about  85  w  m  Hz  ,  U.S.  Standard  Atmosphere  1962, 
and  the  T,  Shimazaki  (1968)  model  which  accounts  for  isomeric  dif¬ 
fusion  due  to  chemical,  temperature,  and  neutral  particle  gradients]. 

The  02  determination  of  the  three  NRL  experiments  agrees  gener¬ 
ally  with  that  of  other  experimenters;  the  latter  data,  however,  do 
not  show  the  rapid  change  of  02  density  with  change  in  altitude  near 
110  km.  This  surprising  feature  of  the  experimental  results  is  not 
reflected  in  either  the  compilations  of  solar  geophysical  data  by 
the  U.S.  Department  of  Commerce  (1966)  or  COSPAR  Standard  Atmospheres 
(CIRA,  1965).  It  may  be  a  temporary  varying  phenomenon. 

Lockey,  Horton  and  Rofe  (1S69)  reported  a  band  of  almost  con¬ 
stant  density  appearing  over  an  interval  of  altitude.  Quessette 
(1970)  reported  data  that  suggest  a  density  step  near  105  km. 
Relatively  few  measurements  of  the  molecular  oxygen  distribution 
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have  been  made  in  the  critical  region  of  the  atmosphere  between  100 
and  120  km.  Consequently,  it  is  possible  that  either  a  constant 
density  over  a  range  of  altitude  (as  described  by  Lockey  et  al., 

1969),  or  a  very  small-scale  height  (rapid  change  of  an  order  of  mag¬ 
nitude  in  short  distance),  as  observed  by  these  experiments,  is  pos¬ 
sible  in  the  terrestrial  atmosphere.  Since  measurements  made  in  this 
region,  in  general,  have  poor  height  definition,  the  profile  over  the 
critical  range  from  100  to  120  km  is  not  clearly  defined. 

3.  Electron  Density  Profiles 

The  electron  density  profiles  shown  in  Fig.  26  indicate  the  fol¬ 
lowing  features:  at  0710  CST  (the  time  at  which  NC7.181  was  launched 
at  Eglin),  the  solar  zenith  angle  was  nearly  the  same  as  that  at  Fo^t 
Churchill  at  1423  CST  (the  time  at  which  NC7.178  was  launched).  How¬ 
ever,  the  E-region  above  120-km  altitude  is  more  intensely  ionized  at 
Fort  Churchill,  presumably  because  of  the  longer  solar  radiation  ex¬ 
posure.  Comparison  with  the  Eglin  data  obtained  at  1357  CST  (time  of 
NC7.180  launch)  reveals  that  at  similar  midday  local  solar  times,  the 
lower  altitude,  smaller  zenii.h-angle  location  (Eglin)  is  characterized 
by  a  greater  degree  of  E- region  ionization. 

4.  Temperature  and  Wind  Profiles 

Temperature  (U.S.  Standard  Atmosphere  1962)  and  wind  profiles 
(Kochanski,  1964,  1966)  are  shown  in  Fig.  29.  Note  that  these  were 
not  measured  at  the  same  time  as  the  absorption  and  ionosonde  measure^ 
ments,  but  represent  the  latitudes  and  reason  in  which  the  NRL  experi¬ 
ments  took  place. 

Under  meteorological  conditions,  the  criterion  for  turbulence 
maintenance  is  commonly  expressed  in  terms  of  the  dimensionless 
Richardson's  Number,  Ri,  which  is  the  ratio  of  the  restoring  buoy¬ 
ancy  force  per  unit  length  of  vertical  displacement,  to  the  square 
of  the  wind  shear. 

Instability  in  the  presence  of  a  density  or  temperature  gradient 
under  gravity  is  implied  by  Richardson  number: 
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Ri  <  Ri 


0.04  -  1 


(12) 


crit 


where  Ricr^t  is  a  critical  value  of  Richardson  number  (Schlichting, 
1960). 

Based  on  Richardson  number  data  (Fig.  29),  it  is  concluded  that 
in  the  lower  thermosphere,  near  the  critical  limit  for  suppression 
by  the  temperature  gradient,  turbulence  may  be  likely,  particularly 
in  the  100-  to  110-km  altitude  range.  The  cessation  of  mixing  due  to 
turbulence  may  determine  the  density  step.  Below  the  step,  tur¬ 
bulent  mixing  facilitates  the  third-body  recombination  of  atomic 
oxygen.  Above  the  step,  the  slower  mixing  by  diffusion  implies  slower 
recombination,  an  enhanced  atomic  oxygen  population,  and  a  reduced 
molecular  oxygen  population. 
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FIGURE  29.  Temperature,  Wind  Speed,  and  Richardson  Numbers 
as  a  Function  of  Height 


61 


5 .  Comparison  of  Homopause  and  Turbopause  Altitudes 


An  altitude  comparison  of  the  step  in  02  density  (called  "homo- 
pause",  signifying  the  altitude  region  below  which  the  atmospheric 
composition  is  essentially  constant)  is  given  in  Table  10.  The  homo¬ 
pause  altitude  is  compared  with  that  of  the  E- layer  bottom  (marked 
by  the  sharp  increase  of  electron  density  observed  by  ionosonde),  and 
with  the  upper  altitude  of  thermosphere  turbulence  (called  "turbo- 
pause")  indicated  by  a  Richardson  number  based  on  representative 
wind  and  temperature  gradients.  The  results  show  agreement  within 
the  uncertainty  of  the  measurement:  the  discrepancy  between  absorp¬ 
tion  and  ionosonde  evidence  of  the  homopause- turbopause  is  about  one 
km  in  the  morning  case  of  the  Eglin  location,  and  almost  exact  for 
«he  midday  case  at  Eglin  and  the  midday  case  at  Fort  Churchill. 

6.  Theory  of  n(0  )  Step  by  Noniscmeric  Diffusion 

Tchen  (1969)  has  given  a  theoretical  explanation  to  account  for 
the  anomalous  reduction  of  0^  density  near  110  km  by  turbulent  dif¬ 
fusion  processes.  He  finds  that  the  density  n(z)  may  vary  as  a  func¬ 
tion  of  altitude  z  as 


n(z)  =  N(z)  f(z)  , 


(13a) 


whttre 


n  T 

N(z)  =  exP 


J-/  dz(Y+H-1) 


(13b) 


and 


f(z)  =  exp 


-  /  -4  *  3?) 


(13c) 
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TABLE  10.  COMPARISON  OF  ALTITUDES  OF  HOMOPAUSE  BY  ABSORPTION  AND  IONOSONDE  OBSERVATIONS 

AND  OF  TURBOPAUSE  BY  REPRESENTATIVE  TEMPERATURES  AND  WINDS 


**Based  on  USSA  1962,  January  temperatures  and  Kochanski  (1964,  1966)  winds,  so  that 
Ricnardson  No.  implies  turbulence. 


(14a) 


also. 


no  =  "<z  =  Zo) 


=  T(Z  =  z0) 


D  -  D  _ 
_  nn  nT 

0  "  T7 


,  +  D  _ 
mol  nT 


In  Eq.  13,  the  effect  of  chemical  change  if,  accounted  for  by 


(14b) 


where  X  =  chemical  reaction  wherein  7  •  X,  is  rate  of  production  of 

a 

density. 

Tchen  finds  two  diffusion  coefficients  effective  in  addition  to 

the  molecular  diffusion  coefficient  D  These  are  D  „  and  D 

mol  nn  nT 

based  on  an  autocorrelation  of  two  density  fluctuations  or  on  a  cross- 
correlation  between  a  density  fluctuation  and  a  temperature  function 
giving  an  eddy  temperature  diffusion  coe  >cient: 


f  dt  n'Mn'M 


(15a) 


and  an  eddy  mass  diffusion  coefficient: 


n'(o)  T'(z) 

- nT - 


(15b) 
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Also,  the  scale  height  is  given  by 


If  the  diffusion  is  "isomeric,"  i.e., 


DnT  =  Dnr>  °r  «  =  °  <16> 

the  solution  (13a)  degenerates  to  (13b)  with  the  reduction  factor 
f (2)  =  1.  The  solution  (13b)  is,  in  such  a  case,  a  smooth  function 
decreasing  with  altitude  in  the  manner  exhibited  by  the  U.S.  Standard 
Atmosphere  (1962). 

In  most  cases  of  turbulent  atmosphere,  the  diffusion  is  not  iso¬ 
meric,  and  DnT  <  Dnn,  so  that  0  >  0. 

Under  such  a  circumstance  of  "nonisomeric"  diffusion,  the  dif¬ 
ference  D  -  D  _  plays  a  role  in  the  distribution  of  n.  One  expects 
nn  ni 

this  effect  to  occur  in  a  region  where  the  turbulent  motion  is  the 
strongest. 

Since  the  reduction  factor  f(z)  does  not  contain  the  chemical 
effects,  it  can  be  calculated  independently.  A  reduction  factor  of 
10  in  an  interval  of  the  scale  height  requires  that  B  =  1.7,  or 
that 


0.35 


nn 


(17) 


a  situation  which  is  not  impossible  in  the  turbulent  atmosphere.  The 
theory  is  elaborated  by  Tchen  (1970). 

The  anomalous  reduction  of  0£  near  110  km  observed  in  the  measured 
data  reported  in  Fig.  28  may  correspond  to  the  effect  of  the  noniso¬ 
meric  diffusion  (Dnn  >  DnT)  by  turbulent  movements. 
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7.  Credibility  of  Geophysical  Determinations 


The  credibility  of  the  geophysical  observations  described  in 
the  previous  paragraphs  is  weakened  by  three  considerations: 

a.  The  step  in  n(02)  at  110  k m  is  not  reported  by  other 
experimenters  and  is  not  traced  reliably  over  its 
whole  range  by  a  single  instrument  in  the  NRL  rocket 
experiments  reported  here. 

b.  Density  values  n(tot)  in  the  range  100  to  120  km 
appear  larger  than  given  by  extrapolation  downward 
of  mass  spectrometer  measurements  above  120  km. 

c.  Temperatures  between  70  and  100  km  appear  to  be 
larger  than  given  by  the  standard  atmosphere  models. 

Considerations  reinforcing  the  credibility  are: 

a.  There  is  redundancy  of  measurements.  Each  wave 
band  channel  in  each  rocket  vehicle  employed  two 
detectors  and  yielded  similar  results.  Three 
rocket  experiments  each  yielded  analogous  results, 
making  allowance  for  expected  geophysical  varia¬ 
tions  due  to  differences  of  time  and  site. 

b.  The  0^  density  data  fits  those  measured  by  other 
experimenters.  At  altitudes  less  than  110  km, 
n(0.,)  is  approximately  22  percent  of  the  "total" 
particle  density,  as  may  be  expected  in  a  homo¬ 
geneous  atmosphere. 

c.  The  base  of  the  E- layer,  as  measured  by  simul¬ 
taneous  ionosonde  at  the  same  time  as  the  absorp¬ 
tion  measurements,  corresponds  in  altitude  with 
the  0^  density  step  (in  two  cases,  closely;  in 
one  case,  approximately) . 

d.  Representative  wind  and  temperature  profiles  indicate 

a  narrow  zone  of  turbulence  may  be  expected  just  below  the 
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02  density  step.  This  local  zone  of  turbulence, 
with  effects  described  analytically  by  Tchen  (1970), 
provides  a  mechanism  for  the  mixing  to  a  homogeneous 
atmosphere  below  the  turbopause-homopause.  Above  the 
turbulence  zone,  where  diffusion  is  the  sole  mixing 
mechanism,  the  atmosphere  becomes  inhomogeneous; 
here,  atomic  oxygen,  produced  by  dissociation  of 
02  but  with  recombination  (principally  by  a  three- 
body  reaction)  reduced,  becomes  a  significant  con¬ 
stituent  of  the  atmosphere.  Below  the  turbulence 
zone,  enhanced  mixing  depopulates  atomic  oxygen  by 
the  more  effective  three-body  recombination. 

e.  Atmospheric  temperatures  are  similar  as  determined 
at  the  same  altitudes  by  different  detectors.  This 
is  evidenced  in  the  slopes  of  density  profiles  for 
both  02  and  "total”  particles,  sine**  the  temperature 
of  the  atmosphere  is,  except  at  the  02  density  step, 
inversely  proportional  to  the  logarithmic  density 
gradient. 

f.  Flux  data  of  measurements  outside  the  effectively 
absorbing  atmosphere  fit  other  observations. 

8.  Summa  ry 

Important  geophysical  features  of  the  experiments  are  summarized 
as  follows: 

a.  A  step  of  02  density,  n(02)(z)  exists  at  altitude 
113  ±  7  km  on  21  February  1966.  There  is  no  com¬ 
parable  slope  discontinuity  of  "total"  particle 
density,  n(tot).  The  step  is  believed  to  mark  the 
homopause. 

b.  The  homopause  is  higher  at  times  of  smaller  solari- 
zation.  Two  experiments  at  like  latitudes,  but 
different  times  of  day,  show  the  homopause  to  be 
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higher  when  the  sun  is  farther  from  the  zenith. 

Two  experiments  at  like  times,  but  at  different 
latitudes,  show  the  homopause  to  be  higher  at  the 
higher  latitude.  Two  experiments  with  comparable 
solar  zenith  angles  show  the  homopause  altitudes 
to  be  similar. 

c.  Homopause  generalizations  derived  by  absorption 
measurements  are  corroborated  by  ionosonde-derived 
electron  density  profiles,  showing  the  marked  slope 
change  identified  as  the  base  of  the  E-layer.  The 
E-layer  base  seems  to  correspond  with  the  homopause 
determined  by  the  step  in  0 2  density,  n(02)(z). 

d.  Wind  and  temperature  profiles,  representative  of  the 
location  and  season  of,  but  not  simultaneously  meas¬ 
ured  with,  the  absorption  and  ionosonde  experiments, 
indicate  that  a  zone  of  strong  wind  shear  tending  to 
turbulence  may  exist  (on  21  February  1966)  just  below 
the  observed  homopause  at  from  101  to  108  km.  The 
upper  bound  of  such  a  shear-induced  turbulence  per¬ 
haps  marks  the  turLopause. 

e.  The  "total”  particle  density  f[n(N2)]  +  N(02)  +  \  n(G)]} 
in  the  100-  to  125-km  altitude  region  exceeds  that  de¬ 
termined  by  mass  spectrometer  measurements,  and  repre¬ 
sented  by  the  U.S.  Standard  Atmosphere  (1962),  by  a 
factor  of  about  two. 

f.  Atmospheric  temperatures,  roughly  determined  to  be 

in  the  150  to  375  °k  range  for  70-  to  150-km  altitudes, 
are  lair  ar  at  altitudes  below  120  km  than  those  of 
the  U.S.  Standard  Atmosphere  (1962)  and  CIRA  (1965). 
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APPENDIX  A 

SOLAR  EMISSION  CONSTANCY 

In  designing  an  experiment  that  spans  an  interval  of  time,  it  is 
desirable  to  ensure,  where  possible,  that  only  intended  variables 
change.  In  the  NC7.178,  180  and  181  experiments,  the  two  variables  of 
the  problem  were  conceived  to  be  the  hour  angle  of  the  sun  and  the 
latitude  of  the  observations.  The  range  of  time  between  the  first 
and  last  launching  was  nearly  seven  hours.  It  was  desirable,  there¬ 
fore,  to  establish  whether  large  transients  other  than  those  due  to 
changing  hour  angle  of  the  sun  and  latitude  of  observation  occurred 
in  the  interval  from  about  0700  to  1400  CST  on  21  February  1966. 
Specifically,  it  was  hoped  that  during  this  interval  of  time  no  sud¬ 
den  change  in  solar  activity  would  occur  to  produce  a  sharp  change  in 
the  normal  quiet  state  of  the  terrestrial  atmosphere. 

Several  types  of  observations  determined  that  solar  emission  was 
comparatively  constant  during  the  period  of  the  experiments.  These 
included  observations  by  satellite,  ionosonde,  and  a  world-wide  riom- 
eter  network. 

1.  SATELLITE  OBSERVATIONS 

The  X^ray  monitoring  satellite  Explorer  30,  which  provided  numer¬ 
ous  measurements  of  solar  X-ray  emissions  during  1966  (Kreplin,  1966), 
has  measured  radiation  in  the  1-8,  8-16,  and  44-60A  wave  ba- ds. 

Data  in  Fig.  A-l  from  McClinton  and  Kreplin  (1967)  show  solar 
flux  measured  by  the  satellite  on  the  day  of  the  NC7.178,  180  and  181 
launches.  This  flux  was  essentially  constant  between  1300  and  2000 
GMT. 
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FIGURE  A-1 .  X-Roy  Flux,  Explorer  30  Satellite,  21  Fobruory  1966 
(McClinton  and  Kreplin,  1967) 


2.  RIOMEIER  OBSERVATIONS 


Riometers  of  a  type  described  by  Little  and  Leinbach  (1959)  are 
located  at  various  stations  throughout  the  world,  approximately  on 
the  75th  meridian  west  and  on  the  geomagnetic  equator.  They  are  oper¬ 
ated  under  the  general  coordination  of  the  APCKL  Laboratory,  Bedford, 
Massachusetts  (Raymond  J.  Cormier,  Principal  Investigator).  The  riom¬ 
eters,  sensitive  to  radiation  of  30  MHz,  indicate  changes  in  radia¬ 
tion  from  sources  outside  the  earth’s  atmosphere,  as  well  as  changes 
in  the  degree  of  ionization  of  the  atmosphere  itself. 

Records  for  a  number  of  the  stations,  plotted  in  local  sidereal 
time  of  21  February  1966,  are  shown  in  Fig.  A- 2. 

When  compared  with  a  plot,  for  each  location,  of  average  signal 
level  for  the  month,  the  data  of  Fig.  A-2  show  that  no  significant 
absorption  events  occurred  on  21  February  1966. 

3.  ICNOSONDE  OBSERVATIONS 

The  degree  of  atmospheric  excitation  may  be  indicated  by  the 
distribution  of  electron  density  in  the  atmosphere.  The  ionosonde 
is  the  most  direct  means  of  measuring  electron  density  of  the  iono¬ 
sphere  from  the  ground. 

Characteristics  of  the  ionosonde  data  taken  hourly  during  Feb. 
ruary  1966,  including  the  day  of  the  launches,  are  described  general¬ 
ly  in  APCRL  Geophysics  and  Space  Data  Bulletin  (1966).  Detailed 
measurements  by  ionosonde  at  the  launch  sites  (Fig.  26)  indicate  a 
normal  diurnal  change  of  the  lower  ionosphere. 

4.  OTHER  OBSERVATIONS  OF  SUN  AND  ATMOSPHERE 

Numerous  other  observations  of  the  sun  and  the  atmosphere  were 
made  (APCRL  Geophysics  Space  Data  Bulletin,  1966,  and  U.S.  Department 
of  Commerce,  ESSA-ITSA  Reports,  1966).  These  imply  the  following 
conclusions  concerning  solar  and  atmospheric  activity  on  21  February 
1966: 
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FIGURE  A-2.  30-MHz  Riometer  Recordings  for  21  February  1966 


No  unusual  disturbances  of  solar  emission  occurred. 

No  strong  solar- indue ed  transients  occurred  in  the 
terrestrial  atmosphere.  This  is  evidenced  by  the 
fact  that  the  APCRL  world  riometer  network  detected 
no  appreciable  30 -MHz  absoi'ption  event. 

A  minor  absorption  event  was  observed  by  a  30 -MHz 
riometer  at  Great  Whale  River,  Canada,  indicating 
the  presence  of  some  auroral  activity  in  the  period 
from  0130  to  2146  UT.  This  interval  includes  the 
times  of  measurements  by  the  NC7.178,  180,  and  181 
experiments. 

Electron-density  profiles  measured  by  ionosonde  at 
the  time  and  launch  location  of  the  NC7.178,  180, 
and  181  experiments  appear  regular. 
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FRACTIONAL  VARIATION  OF  DENSITY  COMPUTATIONS 

The  fractional  variation  or  coefficient  of  variation  (ratio  of 
standard  deviation  of  individual  density  determinations  to  mean 
values)  is  calculated  for  intervals  of  8  km  over  a  part  of  the 
credible  range. 

For  each  altitude  value  within  the  desired  credible  range,  the 
corresponding  density  value  for  one  detector  (Phase  A)  of  some  chan¬ 
nel  is  denoted  by  DA^  and  for  the  other  detector  (Phase  B)  by  DBt  , 
where  i  takes  values  from  1,...,  NP,  and  NP  represents  the  total 
number  of  data  points  within  the  altitude  range  of  the  computation. 
Then: 


AVt  =  (DAt  +  DBt)/2  (Bl) 

For  each  altitude  z^  (t  =1,...,  NP),  N  is  determined  so  that 


'i+N 


£  z.  +  8  £  z 


t+N+1 


and 


79 


ER 


i 


(AV.  -  DA.)2/(N  +  1) 


i+N 

AVj/<N  +  1) 

j  =  t 


(B2) 


A  single  variation  value  is  also  calculated  for  the  entire  cred¬ 
ible  range.  On  this  basis,  density  errors  (ratio  of  standard  deviation 
to  mean  values)  for  experiment  NC7.178  are  shown  in  Fig.  B-l,  for 
NC7.180  in  Fig.  B-2  and  for  NC7.181  in  Fig.  B-3.  Density  variations, 
so  computed,  do  not  include  systematic  errors  in  the  measurement  of 
solar  flux  for  which  estimates  have  been  given  m  the  body  of  the 
report. 
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